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Abstract 

We consider various classes of persistent extragalactic astrophysical sources which 
have been suggested in literature as possible emitters of ultra-high-energy cosmic rays 
(UHECR) . We compare the strength of the claimed correlations by a uniform procedure 
for all classes of sources by making use of the AG AS A, Yakutsk and HiRes stereo data. 
BL Lac type objects correlate with the cosmic rays detected by all three independent 
experiments and are more probably, compared to other astrophysical sources, related 
to the UHECR origin. With the account of the Galactic magnetic field (not possible for 
the HiRes data at the moment), apart of BL Lac type objects, unidentified gamma-ray 
sources may be correlated with AGASA and Yakutsk cosmic rays. 



1 Introduction 

Thousands of cosmic rays with energies higher than 10 19 eV have been detected by various 
experiments. Still, the origin of these energetic particles remains unknown. Deflections of 
charged particles by cosmic magnetic fields are relatively small (though significant) at these 
energies; moreover, presence of a fraction of neutral particles, which propagate rectilinearly, 
is not excluded. This opens a possibility for direct searches of astrophysical sources of the 
high-energy cosmic rays by positional correlations. 

Several difficulties limit the application of this approach. Firstly, at least a large fraction 
of air showers with E > 10 19 eV is believed to be caused by protons, whose trajectories 
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are bend by cosmic magnetic fields by several degrees. Lack of knowledge about these fields 
makes it difficult to account for the deflections while at highest energies, where the deflections 
are believed to be small, low statistics of the cosmic rays prevents one from making definite 
conclusions. 

Secondly, the angular resolution of the cosmic-ray experiments is very poor by the as- 
tronomical standards. The most precise ground array, AGASA, has an average error of 
1.8° in determination of the arrival directions at the highest observed energies, and 2.5° at 
E ~ 10 19 eV p. As a result, an impressive number of astrophysical objects fall in the error 
box of each particular cosmic ray. The situation improved with the availability of the first set 
of the HiRes stereo data with claimed angular resolution of 0.6°. Still, this resolution is 
quite poor for direct identification: dozens of optical, radio and X-ray sources are contained 
in a given circle of 0.6° radius on the Celestial sphere. 

With these complications, statistical methods give the only possible clue to search for 
positional correlations between the cosmic rays and their potential sources. An excess of 
cosmic rays at small angular distances (in particular, at those compared to the experiment's 
angular resolution) from objects of a certain class may suggest that this class of objects 
contains more probably the sources of these cosmic rays than other classes, which do not 
exhibit such an excess. 1 

Many astrophysical sites have been suggested which can host an accelerator powerful 
enough to produce particles with E ~ 10 19 eV and higher. Most of them are extragalactic 
(which is consistent with the lack of global galactic anisotropy of the arrival directions of high- 
energy cosmic rays). Various cosmic accelerators may operate by similar mechanisms, so the 
most efficient way to distinguish between them is based on direct positional correlations. Up 
to now, the correlations with different classes of objects were studied by different methods, 
with different data sets, or were not studied at all. This work aims at filling this gap and 
calculating the chance probabilities of the claimed correlations by a uniform method, which 
allows to compare the strength of correlations and to figure out more probable sources among 
those suggested previously. 

The rest of the paper is organized as follows. In Sec. 2, we discuss the logic of the sta- 
tistical analysis of correlations, describe the methods we use and emphasize several subtle 
points. Sec. 3 lists potential astrophysical sources, briefly recalls motivations for these can- 
didates and previous correlation studies. For each class of sources, we give here references 
for the catalogs we use in the current study. In Sec. 4, we present the results of the cor- 
relation studies for several representative samples of cosmic rays. The conjectures about 
positional correlations were formulated on the basis of the cosmic rays observed by AGASA, 
and sometimes by the Yakutsk experiment. We use the samples of cosmic rays with energies 
4 • 10 19 < E < 10 20 eV observed by AGASA (49 rays) and by AGASA and Yakutsk (59 rays) 

lr The background of random coincidences is always present, and one cannot claim that one particular 
object is a source of a given cosmic ray based on these statistical arguments only. 
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to treat the suggested sources uniformly and to reformulate the conjectures. We consider 
independently neutral and proton primaries and in the latter case, we account for proton 
deflection in the Galactic magnetic field. Then, we test the same conjectures with recently 
published HiRes stereoscopic data (271 cosmic ray with E > 10 19 eV), assuming only neutral 
primaries, since energies of the HiRes events are unpublished. Finally, we look for correla- 
tions of suggested sources with the cosmic rays of the highest (E > 10 20 eV) energies. Sec. 5 
contains conclusions and discussion. 

2 The correlation analysis 

The general logic of the correlation studies by statistical methods is to test the hypothesis 
that the distribution of the observed arrival directions of the cosmic rays is isotropic, and that 
there is no excess of the cosmic rays at small angular distances 2 from the sources listed in a 
given catalog (and not to test the hypothesis that the sources from the list emit cosmic rays). 
To this end, a large number of cosmic-ray events are simulated which arrive at the Earth 
isotropically. Some part of them is "accepted", following the same rules which determine the 
acceptance of a given cosmic-ray experiment. The angular distances between the simulated 
events and the astrophysical objects from the catalog are calculated and their distribution is 
compared to the distribution of distances between the real cosmic rays and the same objects. 
The less the probability that the real distribution is a statistical fluctuation of one calculated 
for isotropic cosmic rays, the stronger the indication that the objects from that particular 
list may be the cosmic-ray sources. 

Another way of reasoning might be to simulate the positions of the sources and to compare 
them with real cosmic rays, or to simulate both. However, the distribution of the objects 
of a particular class in the Universe follows specific laws and is not easy to be simulated. 
The selection biases of discovery, observation, and inclusion in catalogs add more to this 
problem. On the other hand, the directional acceptance of the cosmic-ray detectors is quite 
well understood and the distribution of the observed cosmic rays agrees with the theoretical 
acceptance reasonably well. That is why we simulate the cosmic rays and not the sources. 

To compare the real and simulated distributions of angular distances, we follow a conven- 
tional procedure used in correlation studies. One takes a catalog of n r cosmic rays observed 
by a given experiment and calculates the number v{$) of events falling within the angular 
distance S from the suggested sources. The value of 5 should be of order of the experiment's 
angular resolution. The actual value at which the real effect is better seen may be deter- 
mined by simulations described in Ref. [3] for AGASA and Yakutsk and in Ref. for HiRes 
stereo. In all cases, the simulation gives the value of 5 somewhat larger than the angular 
resolution. The reason for this smoothing is quite clear: a significant part of the observed 
cosmic rays is not related to known sources included in catalogs. The Monte Carlo code 
2 That is, at those of order of the experimental angular resolution. 
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generates then a random direction in the visible half of the Celestial sphere, determined by 
the azimuth angle A and the zenith angle Z, as well as the sidereal time of the event arrival, 
s. Then the code either accepts (with the probability P(A, Z, s)) or rejects the event. The 
procedure is repeated until iV sets, each of n r accepted events, are generated, where N is 
large: we use N = 10 4 sets in our study except for the cases of a very strong signal when we 
use N = 10 5 . For the ground arrays, the probability P(A, Z, s) = cos Z\ for the HiRes exper- 
iment in stereo mode, P(A, Z, s) = Pa(A)P z (Z)P s (s) and the functions Pa,Pz,P s are given 
in Ref. [Sj. For each mock set i = 1, . . . N, the quantity v%{8) is calculated. The probability 
that the observed distribution of cosmic rays with respect to the sources is a fluctuation of 
a random distribution is estimated (for each particular pair of samples of cosmic rays and 
astrophysical objects) simply by 

number of mock sets withz/j(<5) > v{$) 
= N ' 

It is worth noting that the value of V should not be, in general, interpreted as the statistically 
correct probability if subsequent tries with different catalogs were made (see Sec. 14. 3|) and if 
the purpose of a study is to test the null-hypothesis that there are no positional correlations 
between cosmic rays and any of the astronomical catalogs considered. 

For the cosmic-ray experiments with relatively good angular resolution, and assuming 
charged primaries, the signal may be smeared out by deflections in the magnetic fields. One 
may attempt to correct the arrival directions for bending in the Galactic magnetic field 
(GMF), assuming a particular GMF model. In this study, when applicable, we assume the 
GMF model used in Ref. jHj and proton-like primaries (electric charge +1). The knowledge of 
particle energies is crucial for the correction for GMF and 30% errors in energy determination 
may result in additional increasing of the value of 5. All simulated events are subject to the 
same correction as the real ones; to this end, the MC code assigns randomly the energies of 
n r observed events to n r simulated events in each of A" mock sets. 

The cuts imposed by hand on the samples of cosmic rays and possible sources may 
result in essential lowering of the probability V. One may impose these cuts a posteriori, to 
maximize the signal, and to compensate the probability by a penalty factor (see, for instance, 
Ref. jS]). In this study, we do not impose any additional cuts on the catalogs taken from 
original works. All cuts in astronomical catalogs (e.g. the z < 0.01 cut in a sample of Seyfert 
galaxies or several cuts in a sample of dead QSO candidates, see Sec. 3) were made in the 
original works where these candidate objects had been suggested. We refer the reader to 
these works for motivation of the cuts. In the catalogs of cosmic rays, the cuts in energies and 
zenith angles were made a priory as well and were determined by the availability of data. We 
note here that making cuts which maximize the signal is very useful for further detalization 
of the class of cosmic-ray emitters but postpone this subject to subsequent studies. 

One subtle point in this kind of a study concerns the value of the bin size 5. The 
probability V depends strongly on 5 (see plots in Sec. 4). Hence, choice of 5 may be considered 
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as a cut which requires its own penalty factor. However, varying 5 at scales smaller than 
the angular resolution of a detector, one can easily obtain strong unphysical correlations: 
one single coincidence within 0.1° of a source and a cosmic ray observed by a detector with 
the angular resolution of 2° results in a drop of V of order (0.1/2) 2 = 2.5 • 10~ 3 compared 
to one coincidence within 2°, while physically the two cases have the same meaning (see 
Fig. 3 in Sec. 4.4). Account of the correlations of this kind in simulated sets may result 
in strong overestimation of the penalty factor if one chooses to determine 5 a posteriori to 
maximize the signal and then to compensate for this choice by introducing penalties. In this 
study, we fix 5 a priori following the simulations of Refs. E] • These values are of order 
of the experiment's angular resolution: ps 2.5° for AGASA (2.8° for AGASA with GMF 
correction) 3 and ~ 0.8° for HiRes stereo. We will see that in the cases when a strong signal 
is seen, the minimal V often corresponds to a different but close value of S. 

Since we are interested in a comparative study of various astrophysical sources only, 
hence in relative values of probabilities, and not in testing a hypothesis that the objects of 
a particular class emit cosmic rays, this approach is sufficient for our purposes. However, a 
more adequate method to obtain absolute values of probabilities should avoid the choice of 5 
at all. Indeed, there exist statistical methods to compare distributions without use of explicit 
binning. Probably the best known one is the Kolmogorov-Smirnov (KS) test; however, in our 
case it is not applicable for two reasons. Firstly, the KS test does not account for the value of 
experimental errors. Secondly, its sensitivity to differences between distributions close to the 
end points of the interval at which the distributions are determined is very poor, while the 
positional correlations appear at values of the angular distances compared to the detector's 
resolution and much less than the maximal possible angular distance of 180°, that is very 
close to the zero end point of the interval [0°, 180°]. One may attempt to cut the interval at 
some angle A in order to shift the effect towards the center of the interval; we checked that 
the result is very sensitive to A in this case. Methods other than the KS test exist which 
are more sensitive at the ends of an interval, but, as, for instance, the Anderson-Darling 
method, they are directly implemented only for the simplest distributions. As a result, this 
approach awaits further studies. 

3 Potential sources 

Various classes of potential astrophysical accelerators are discussed, for instance, in Ref. [Zj. 
We tried to include in this study all sources mentioned in the literature as candidate sources 
of ultra-high-energy cosmic rays, subject to two criteria: we require the objects to be per- 
sistent (gamma-ray bursts are excluded in this way) since the comparison of significance 
can hardly be performed between persistent sources and flares, and we consider only extra- 
galactic objects since no galactic anisotropy is seen with the impressive statistics of events 
3 We use the same values for the AGASA-dominated joint samples of cosmic rays. 
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at E > 10 19 eV. For the case of unidentified gamma-ray sources, to safely exclude Galactic 
populations, we require the galactic latitude |fe| > 10°. For comparison, we present also 
results without this cut. 

3.1 Flat spectrum radio quasars 

Powerful active galactic nuclei, where particles can be accelerated up to extreme energies, are 
often considered as the best astrophysically motivated sources of energetic cosmic rays. The 
significant correlations (V = 5 ■ 10 -3 ) of the highest-energy (at least one standard deviation 
above 8 • 10 19 eV) events observed by Fly's Eye, AGASA and Haverah Park and published 
before 1998 with flat-spectrum radio-loud quasars (FSRQs) were first claimed in Ref. jH] 
and considered as an argument in favour of GZK-violating new composite particles. Since 
then, these correlations were reanalyzed in various papers [01 EI| with different results. The 
original sample consisted of quasars from the Bonn catalog ^Tj of radio sources with 5 GHz 
flux exceeding 1 Jy and radio spectral index a > —0.5. We considered both the original 
sample and the sample selected by the same criteria from the new 11th edition of the Veron 
catalog of quasars [T2j . 

Catalog 1: Kuehr FSRQs 

Original catalog: Kuehr et. al., 1981 jTTj. 
Cuts: class=QSO, spectral index a > —0.5. 
Number of objects: 201. 

Catalog 2: Veron FSRQs 

Original catalog: Veron-Cetty and Veron 2003 ^21, Table 1 (quasars). 

Cuts: 6 cm flux F§ > 1 Jy, spectral index a > —0.5 (calculated from fluxes at 6 and 
11 cm). 

Number of objects: 211. 

3.2 Seyfert galaxies 

This subclass of active galaxies was considered in Refs. ^3] where physical motivations were 
discussed and correlations (with the probability of chance coincidence as low as 2 • 10~ 5 ) 
with AGASA and Yakutsk cosmic rays with energies E > 4 • 10 19 were claimed for several 
samples with fine-tuned cuts in the galactic latitude. In Ref. the catalog ^1] and the 
10th edition of Veron catalog [TH] were used and the distance cut z < 0.01 was imposed. We 
use both the catalog [Hj and the updated Veron catalog [12 j with the same cut. 

Catalog 3: Lipovetsky Seyferts 

Original catalog: Lipovetsky et. al. |14j . 
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Cuts: z < 0.01. 
Number of objects: 136 

Catalog 4: Veron Seyferts 

Original catalog: Veron-Cetty and Veron 2003 J2]> Table 3 (active galaxies) 
Cuts: classification starts with S (all Seyferts), z < 0.01. 
Number of objects: 190. 

3.3 BL Lac type objects 

Active galaxies often have jets fueled by the nuclei and consisting of relativistic plasma fired 
in strongly collimated cones. When one of the jets is directed towards the observer, this 
galaxy is seen as a blazar. These objects are naturally considered as good candidates for the 
sources of the cosmic rays which reached the Earth. BL Lac type objects constitute a subclass 
of blazars distinguished by their spectral properties. Different authors use different criteria 
for this distinction. In particular, one may either require the absence of the emission lines or 
put a constraint on their width. These criteria are discussed, for instance, in Ref. fB] and 
are implemented in the catalog ^2] where Table 2 contains both confirmed BL Lac's in the 
strict sense of the word and possible BL Lac's, notably those classified as high-polarization 
(HP) blazars. 

The first claim for correlation of cosmic rays with confirmed BL Lac's (V ~ 6 • 10~ 5 ) was 
made |3] for the maximally autocorrelated sample of cosmic rays detected by the AGASA 
and Yakutsk experiments. The sample of BL Lacs was chosen from the 9th edition of the 
Veron catalog [TJj by the cuts on the visual magnitude, the radio flux and the redshift. The 
sample was chosen in such a way that the signal was maximized; cited probability includes 
the penalty factor for cut adjustment j^j. 

As discussed below (see Sec. 3.8), gamma-ray emission should accompany emission and 
propagation of UHECR. It was observed indeed that the BLL which correlate with UHECR 
in the resulting subsample of Ref. [3j are classified as EGRET sources in astronomical 
databases. Imposing the single criterion - positional coincidences of objects in Veron 
and EGRET [22] catalogues - leaves 14 objects and this subsample also exhibits very low 
formal values of V ~ 3 • 10~ 7 with the maximally autocorrelated set of UHECR ^H]- R is 
difficult to calculate the significance here since the set of 14 BLL overlaps with previous set of 
22 BLL of Ref. jHj . For our comparative study we have chosen one of these subsamples (the 
other one shows very similar values of V in all tests we carry). Namely, we have chosen single 
physical criterion of BLL being the EGRET source for our catalog selection. This condition 
is more transparent and is easier to implement with updated most recent catalogue of AGN 
|12j . which we use uniformly throughout our comparative study. 

A much more general sample of objects — all confirmed BL Lac's brighter than 18 m from 
the 10th edition ^3] of the Veron catalog — was used, with a positive result, for the searches 
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for correlations with all published AG AS A events (see Refs. 0113 1201)- in particular, it has 
been noted |H] that the account of the GMF improves the correlations. 

Very recently, the correlations of all three previously selected samples of the BL Lac 
type objects were tested |lj with the HiRes stereo data with the positive result - probability 
V — 4 • 1CT 4 - for the sample [HI HH1 12D| of bright (m < 18) confirmed BL Lac's. 

The correlations with BL Lacs were reported also in Ref. [2~Tj . 

For this study, we use the updated catalog [12] and do not impose any new cuts. We 
consider separately "strict" and high-polarization BL Lac's. We also consider the subsamples 
of the new catalog [T2j selected by the criteria of Refs. jB] and [T%] . The results for the 
old ^3 EI and new J2] catalogs are similar. 

Catalog 5: Confirmed Veron BL Lac's 

Original catalog: Veron-Cetty and Veron 2003 [T2"] . Table 2 (BL Lac type objects). 
Cuts: classification=BL (confirmed BL Lac's). 
Number of objects: 491. 

Catalog 6: HP Veron BL Lac's 

Original catalog: Veron-Cetty and Veron 2003 [12], Table 2 (BL Lac type objects). 
Cuts: classification=HP (high-polarization BL Lac's) 
Number of objects: 68. 

Catalog 7: Bright confirmed Veron BL Lac's 

Original catalog: Veron-Cetty and Veron 2003 [12, Table 2 (BL Lac type objects). 
Cuts: classification=BL (confirmed BL Lac's), visual magnitude m < 18. 
Number of objects: 178. 

Catalog 8: Possible EGRET BL Lac's 

Original catalog: Veron-Cetty and Veron 2003 [12], Table 2 (BL Lac type objects). 

Cuts: classification=BL (confirmed BL Lac), angular distance to the nearest source from 
the 3EG catalog [22] less than 2rgs, where rgs is the angular radius of the circle containing 
the same solid angle as the 95% confidence contour of that EGRET source. 

Number of objects: 30. 

3.4 Radio galaxies 

Often intrinsically as luminous as blazars, similar objects seen at large angles to the jet 
appear much fainter and are called radio galaxies. Their classification has been introduced 
in Ref. 23J; the FR I type galaxies are thought to be the counterparts of the BL Lac's 
while more powerful FR II type objects may be the counterparts of other blazars. The lobes 
(clouds of matter fueled by the jets) and the hot spots (regions where the jet hits the lobe) 
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of FR II type galaxies are efficient acceleration sites [21]. According to simple estimates 
(see, for instance, Ref. [2H| ), the hot spots are the only known places in the Universe where 
protons can be accelerated to energies higher than 10 20 eV by shock acceleration. The shocks 
move in various directions, the magnetic fields in the source are very strong and, as a result, 
the accelerated particles may escape the accelerator not necessarily along the jet. We are 
not aware of any previous correlation searches for the radio galaxies. Here, we use a catalog 
of radio galaxies with their FR classifications published in Ref. [2E] • 

Catalog 9: FR I radio galaxies 

Original catalog: Zirbel and Baum [2E], Table 11. 
Cuts: classification^ or Ig (confirmed FR I). 
Number of objects: 73. 

Catalog 10: FR II radio galaxies 

Original catalog: Zirbel and Baum [2B], Table 11. 
Cuts: classification=II, Ilg, Iln, or lid (confirmed FR II). 
Number of objects: 116. 

3.5 Colliding galaxies 

Regions of strong shocks are obviously present in merging galaxy systems. Based on this 
fact and on several positional correlations (without a detailed statistical analysis), nearby 
colliding galaxies were suggested as possible sources of the highest-energy cosmic rays (see, 
for instance, Refs. PUI2I])- Our sample consists of objects classified as "pair in contact" in 
the latest update of the Vorontsov-Velyaminov catalog of interacting galaxies |28j . 

Catalog 11: Colliding galaxies 

Original catalog: Vorontsov-Velyaminov et. al. [28.. 
Cuts: classification starts with PK (pair in contact). 
Number of objects: 454. 

3.6 Starburst and luminous infrared galaxies 

In a region of large-scale star formation (a starburst region), a strong shock wave may 
develop which provides the conditions required for the efficient particle acceleration. Two- 
step acceleration process in these regions may result [2H| in nuclei energies of 10 20 eV and 
higher. When two colliding galaxies finally merge, they may produce an active starburst 
region near the common center; this appears to the observer as a luminous infrared galaxy 
(LIG). The possibility that the latter can be a source of highest-energy cosmic rays has been 
discussed in Ref. j3U|- Given a large number of sources, the correlation study was possible 
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only in the sense that the large-scale distribution of cosmic rays and LIGs are correlated; 
these studies did not give definitive results [3U1 |3T]. In Ref. |7j, the most probable candidates 
among LIGs were pointed out: those listed in the PDS [3j2; and HCN jHH| catalogs. It has been 
pointed out in Ref. [7J that it is worth considering the list of Ref. (HI] in correlation studies. 
The latter list consists of those LIGs for which the expected gamma-ray flux (calculated in 
Ref. [31]) exceeds the sensitivity limit of GLAST, a future gamma-ray telescope. Though 
it is not perfectly clear why this particular GLAST-related cut is important for ultra-high- 
energy cosmic ray production, we include the sample of Ref. [HI] in our study together with 
the PDS and HCN samples. 

Catalog 12: PDS starbursts 

Original catalog: The Pico dos Dias Survey [32], Table 1. 
Cuts: no. 

Number of objects: 203. 

Catalog 13: HCN luminous infrared galaxies 
Original catalog: HCN survey [SB], Table 1. 
Cuts: no. 

Number of objects: 53. 

Catalog 14: Selected luminous galaxies 

Original catalog: HCN [33] and PDS [32] surveys, compiled in Ref. [31]. 

Cuts: calculated in Ref. jS] expected gamma-ray flux exceeds the GLAST sensitivity. 

Number of objects: 16. 

3.7 Dead quasar candidates in luminous nearby galaxies. 

The observation of the cosmic rays with extremely high energies, at which the propagation 
length of a particle is limited due to its interaction with the cosmic background radiation, 
suggested that their sources should be located in the neighborhood of our Galaxy. Since no 
objects, powerful enough, are explicitly seen at these distances, it has been proposed that 
the highest-energy cosmic rays may be accelerated in quasar remnants which do not emit 
strong radiation observed by astronomical, in particular radio, telescopes [SSI EE]. A typical 
feature of a source of this kind is an extremely heavy black hole in the nuclei of an elliptic 
galaxy. In Ref. [35j, a list of possible candidate sources has been suggested; we included it in 
our study. According to a refined study of Ref. [37], a typical dead quasar candidate which 
has chances to correlate with the observed cosmic rays should satisfy the following criteria: 

1) absolute magnitude M a b s < — 21; 

2) redshift z < 0.01; 

3) morphological type t < — 3 (elliptical); 
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4) galactic latitude |6| > 20°; 

5) the object belongs to a group with richness not exceeding 50. 

The authors of Ref. jHZj applied these five cuts to the Nearby Optical Galaxy catalog [38J. 
The resulting sample consisted of 12 objects. We did not succeed in reproducing their catalog 
of 12 objects by the same procedure, probably due to the fact that the cut (5) is quite vague 
and Ref. [37] used unpublished data for its implementation. In our study, we skipped this 
cut and left only first four. We applied them to the LEDA database |HH] which is more 
complete than Ref. [SB]- Still, not all 12 objects of Ref. jHJJ entered the resulting sample. 
Taking the original 12-object sample, we reproduced the results for the AGASA experiment 
(without the GMF correction) presented in Ref. [31] • This sample exhibits no correlations 
either upon account of the GMF or with the HiRes data. 

Catalog 15: dead QSO candidates from Ref. [35] 

Original catalog: sample of nearby galaxies with massive dark objects, Ref. [IH] . 
Cuts: central object mass M > 1O 9 M . 
Number of objects: 14. 

Catalog 16: dead QSO candidates with cuts of Ref. |37j 
Original catalog: LEDA database [3~§j . 
Cuts: M abs < -21, z < 0.01, t < -3, \b\ > 20°. 
Number of objects: 23. 

3.8 Gamma ray sources 

Both acceleration of particles and their propagation through the cosmic background radiation 
are accompanied by emission of energetic photons. These photons, in turn, lose energy in 
interactions with background photons (Hj and arrive at the Earth with energies of order 
100 MeV - 100 GeV, that is in the range observed by the EGRET instrument. Hence, 
the gamma-ray emission is very probable to accompany the cosmic ray emission (see, for 
instance, Refs. [TBI E21 HH] ) ■ The TeV emission has been suggested as a characteristic feature 
of the AGNs which accelerate protons up to highest energies but are quiet in other energy 
ranges, analogues of the "dead quasars" [4*4] . 

Correlations of the blazars present among the identified sources in the Third EGRET 
catalog [22] were studied in Ref. JU] with a negative result. Most of the EGRET blazars are 
strong sources of the FSRQ or HP type (see Sec. 3.1, 3.3). In Ref. [IB], strong correlations 
(V ~ 10 -4 ) were found between the sample of AGASA and Yakutsk cosmic rays and uniden- 
tified EGRET sources. An important ingredient for these correlations was the inclusion of 
the GMF correction. Another correlating sample of EGRET sources studied in Ref. ^B] 
consisted of 14 sources which have a confirmed BL Lac nearby (see Catalog 8 in Sec. 3.3). 
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Here, we separately study the samples of EGRET sources (i) identified as confirmed 
AGNs in the 3EG catalog and (ii) unidentified (outside the Galactic plane). We also include 
a sample of GeV sources laying outside the Galactic plane and a small sample of extragalactic 
TeV sources. For comparison, we present the results for unidentified EGRET and GeV 
sources without the galactic latitude cut. 

Catalog 17: EGRET blazars 

Original catalog: The 3rd EGRET catalog [22]- 
Cuts: identification=A (confirmed AGN). 
Number of objects: 67. 

Catalog 18: Unidentified EGRET sources outside the Galactic plane 
Original catalog: The 3rd EGRET catalog (221 - 
Cuts: identification void, galactic latitude |6| > 10°. 
Number of objects: 96 

Catalog 18a: Unidentified EGRET sources 
Original catalog: The 3rd EGRET catalog [22]- 
Cuts: identification void. 
Number of objects: 170 

Catalog 19: GeV sources outside the Galactic plane 
Original catalog: Lamb, Macomb 
Cuts: galactic latitude \b\ > 10°. 
Number of objects: 26 

Catalog 19a: GeV sources 

Original catalog: Lamb, Macomb 
Cuts: no. 

Number of objects: 52. 

Catalog 20: Extragalactic TeV sources 

Original catalog: List of confirmed TeV sources presented in Ref. [46J. 
Cuts: identification with extragalactic objects 4 . 
Number of objects: 10. 

4 The only TeV source present in the list of Ref. 0^] and lacking firm identification is most probably 
Galactic and is not included. 
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H 


HiRes stereo 




> 10 19 eV 


m 


271 
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20+ 


All Northern 




> 10 20 eV 


See Ref. gHj 
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Table 1: The cosmic-ray samples. Notes: 1) The updated sample of events observed by the 
Yakutsk array before 1995 with zenith angles less than 45°. 2) The arrival directions were 
published in the form of a PS plot in Hammer projection; the energies of the events are 
unpublished, so the correction for the GMF is not possible. 3) The sample is described in 
detail in Ref. |48j . 

4 Results 

For the catalogs described in Sec. 3, we estimated, as described in Sec. 2, the probability that 
the observed number of cosmic-ray events close to the astrophysical objects is a fluctuation 
of the random distribution. 

4.1 Cosmic-ray samples 

The samples of cosmic rays we use are determined mostly by the availability of data. Most 
of the previous studies were based on AGASA and (sometimes) Yakutsk events. The HiRes 
stereo sample became recently available. Motivated by the fact that several experimental 
groups published only the data on the events with energies higher than 10 20 eV and keeping in 
mind the possibility that the "sub-GZK" and "super-GZK" events may originate in different 
classes of objects, we consider separately the events with E > 10 20 eV. We describe the 
cosmic-ray sets in Table 1. 

4.2 Results of the comparative study 

The values of the probability V calculated following Sec. 2 are presented in Table 2. One 
can make following conclusions out of these Tables. 

1. At present, cosmic rays with E > 10 20 eV do not correlate with any astrophysical 
objects. This does not exclude, however, any real connection because of low statistics. 

2. The strongest correlations seen in AGASA and Yakutsk datasets at energies below 
10 20 eV are exhibited by the catalogs 3, 5, 7, 8 and 18. 
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Table 2: Probability P calculated for the catalogs of astrophysical objects described in Sec. 3 
and the samples of cosmic rays described in Sec. 4.1. For convenience, with the number of 
a catalog, a reference to a section of this paper is given as follows: FSRQ - Sec. 3.1, Sy - 
Seyfert galaxies - Sec. 3.2, BL - BL Lac's - Sec. 3.3, RG - radio galaxies - Sec. 3.4, coll. - 
colliding galaxies - Sec. 3.5, LIG - starburst and luminous IR galaxies - Sec. 3.6, dQSO - 
dead QSO candidates - Sec. 3.7, 7 - gamma-ray sources - Sec. 3.8. The columns with the 
results obtained with correction of the arrival directions for GMF are marked by asterisks. 
The values of V < 4.6- 10 -2 (which would correspond to a 2a effect for the Gaussian statistics) 
are given in bold face. The values V < 2.7 • 10~ 3 (would-be 3<r) are underlined. Each value 
of V has been calculated without accounting for other tries and thus can be interpreted as 
a relative probability with respect to other catalogues. 
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(a) The effect in Catalog 3 (Lipovetsky Seyferts) is dominated by 5 sources close to 
the direction of the AGASA triplet. These sources contribute 8 out of 13 observed 
coincidences with AGASA cosmic rays within 2.8°. The effect is not present for a 
sample of Seyfert galaxies from the Veron catalog. The correlations are observed 
with the account of the correction for GMF only. 

(b) The correlations of Catalog 5 (all confirmed BL Lac's) and Catalog 7 (bright 
confirmed BL Lac's) with AGASA and Yakutsk cosmic rays are improved with 
the correction for GMF. 

(c) The Catalog 8 represents an update of the sample of Ref. ^H] of BL Lac's which 
are possibly associated with EGRET sources. It contains 30 objects, compared 
to 14 original. The values of probabilities presented here are higher compared to 
those given in Ref. • The difference is due to the fact that for the comparative 
study, we use here different sets of cosmic rays. With the update 30-objects 
catalog and the same maximally autocorrelated sample of cosmic rays as used 
in Ref. ^Bj, we obtain the results very similar to those presented there for the 
14-object sample. 

(d) Strong correlations of cosmic rays observed by the ground arrays are seen with 
unidentified EGRET sources (catalog 18) outside the Galactic plane. Generally, 
correlated sources were detected neither in GeV nor in TeV. The physical prop- 
erties of the correlated sources and their identification require further study. The 
correlations with unidentified EGRET sources persist even when a large Galactic 
population is included (Catalog 18a). Note in passing that on the base of the 
correlations with 14 BL Lac objects associated with EGRET sources it was sug- 
gested PH] that at least some of the unidentified EGRET sources, correlated with 
UHECR, are BL Lac type objects still unobserved in other wave bands. 

3. Results on AGASA and ( AG ASA+ Yakutsk) samples are consistent with each other. 

4. The corrections for the deflection in GMF are important. A full test of the correlations 
found in the AGASA and Yakutsk data with the new HiRes stereo data will be possible 
only after publication of the energies of HiRes events. One should note however that 
the poor knowledge of the GMF results in errors in corrected arrival directions which 
exceed the angular resolution of the HiRes experiment in the stereo mode, so the effect, 
if present, may be smoothed considerably. 

5. Only BL Lac type objects (catalogs 5, 7) exhibit correlations both with the HiRes stereo 
data set and the AGASA-dominated samples. This fact is very interesting and deserves 
further discussion; see also Sec. 14.31 To the correlations of HiRes cosmic rays with the 
TeV sources (catalog 20), two objects contribute, and both are BL Lacs. Other catalogs 
are not, however, rejected. For instance, the correlations with unidentified gamma-ray 



15 



sources are mostly seen for the proton primaries when the deflection of arrival directions 
in GMF were taken into account; a similar analysis for the HiRes sample is impossible 
at the moment because the energies of the cosmic rays are unpublished. Moreover, the 
angular resolution of HiRes stereo is better than uncertainties both in the positions 
of the EGRET sources (which means that the ultra-high-energy astronomy is actually 
possible) and in the corrections for GMF. 

4.3 Significance of multiple tries 

When the purpose of a study is to estimate the probability that some particular class of 
objects is related to a sort of events (UHECRs in our case), multiple tries of various catalogs 
are not welcome. Just by chance, one expects the probability of order 1/M to be found 
among M independent catalogs of objects tested for correlations with a single cosmic-ray 
sample. Thus, a factor of M is usually required as a penalty in this case. When the tries 
are not independent, the penalty factor is more difficult to estimate. The purpose of our 
study is, however, purely comparative: we test the conjectures made in previous works by a 
uniform method. To compare the probabilities, we thus do not care about any general factor 
by which all the numbers should be multiplied, and this penalty factor in particular. 

On the other hand, one notices that quite low V are found for the Catalog 7 with the 
cosmic rays observed by different experiments. While this is not our main goal, one may 
wonder how probable is to obtain V < 7 ■ 10~ 4 in independent experiments just by chance. In 
this way, one would test the hypothesis that there are no catalogs correlating with UHECR 
samples among all we list and all the correlations seen are nothing but statistical fluctuations 
due to multiple tries. We performed a simple estimative simulation to this end. We took 
the declinations of objects present in our 21 catalogs and randomized the right ascensions 
many times. In this way, we obtained 32000 sets, each of 21 catalog with the same numbers 
of objects and the same distributions in declination as the real ones. All 21 • 32000 ~ 7 • 10 5 
simulated catalogs were subject to the same procedure as the real catalogs. Only in one 
case, a simulated catalog exhibited V < 7 ■ 10 -4 for the HiRes sample and one of four other 
(interdependent) samples simultaneously. We thus estimate the probability that the results 
for the Catalog 7 are a result of a chance fluctuation as 1/32000 ~ 3 • 10~ 5 . 

4.4 Dependence on 5 

As it has been already noted, the values of the bin sizes we fixed for the calculation of V 
were determined by simulations. In some cases, lower values of V were obtained for other 
choices of 5. In figures 1 and 2, we present the 5 dependence of V for two interesting cases: 
unidentified EGRET sources with AGASA and Yakutsk cosmic rays and BL Lac's with 
HiRes cosmic rays. For comparison, we present on the same plots the V(S) dependence for 
a typical catalog which does not exhibit any correlation. An example of the effect described 



16 




01 23456789 



Figure 1: V{$) for catalog 18 (unidentified EGRET sources outside the Galactic plane, lower 
(red) curve) and catalog 4 (nearby Veron Seyferts, upper (green) curve), cosmic-ray sample 
A+Y (AGASA and Yakutsk), correction for GMF included. 




2 4 6 8 



Figure 2: V(S) for catalog 7 (bright confirmed BL Lac's, lower (red) curve) and catalog 4 
(nearby Veron Seyferts, upper (green) curve), cosmic-ray sample H (HiRes stereo). 
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Figure 3: Curve with a shallow minimum (green): V(5) for catalog 9 (colliding galaxies), 
sample A, GMF included. The dip at 5 = 0.3° is probably unphysical and diminishes for 
5 = 2.8° expected from the experiment angular resolution. Curve with deeper minimum 
(red): catalog 20 (extragalactic TeV sources), sample H. For this sample, V(0.2°) < 10~ 4 
while physically meaningful P(0.8°) ~ 10~ 2 . 

in Sec. 2 - a random coincidence at very small angles - is presented at Fig. 3. 

5 Conclusions 

We performed a comparative study of correlations between the arrival directions of the cosmic 
rays observed by AGASA, Yakutsk and HiRes experiments with 21 catalogs of astronomical 
objects suggested by different authors as possible UHECR sources. No a posteriori cuts were 
imposed on the catalogs. The results are described in Sec. 4.2 but should be interpreted 
with great care because the angular resolution of the cosmic-ray experiments is poor, cosmic 
magnetic fields are known not quite well and the statistics is not overwhelming. 

The analysis of the AGASA and Yakutsk samples demonstrates the presence of cor- 
relations with unidentified EGRET sources, BL Lac's and nearby Seyfert galaxies. The 
corrections for the Galactic magnetic field are very important for these effects. 

Only BL Lac type objects exhibit correlations with the HiRes stereo sample as well. 
With the uncut catalog of 491 confirmed BL Lac's and the set of 271 cosmic rays with 
energies E > 10 19 eV, we obtained the probability of 0.04 for these correlations to result 
from a random coincidence. For the updated catalog of 178 bright (m < 18) BL Lac's, 
the probability is as low as 5 • 10~ 4 . For the AGASA sample (with the GMF correction), 
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V = 4 • 10~ 4 for the same catalog; for AGASA and Yakutsk V = 7 • 10~ 4 . The probability 
that one of 21 catalogs exhibits, by chance, V < 7 ■ 10~ 4 both for the HiRes sample and one 
of four other samples is about 3 ■ 1CT 5 . The correlations are present at the angular distances 
of order HiRes angular resolution and much less than the expected deflections of charged 
particles in cosmic magnetic fields. This suggests that a fraction of the cosmic-ray flux at 
energies E > 10 19 eV is carried by neutral primaries 0]. 

Clearly, more detailed information about the cosmic rays observed by HiRes would be 
very important for understanding the effects discussed here. Knowledge of the energies of 
particles would allow to reconstruct (to some accuracy) their deflections in the Galactic 
magnetic field. A study of shower profiles might help to understand what kind of neutral 
particles may be present in the cosmic rays and contribute to correlations. Finally, if, 
for example, a neutral primary is a photon emitted by a blazar, the arrival time of the 
corresponding air shower may coincide with a temporary change in the state of the source. 
This would give much stronger support to conjectures about particular cosmic-ray sources 
than the statistical studies of the arrival directions. 

Further detailed study will specify which particular physically distinguished subclasses of 
BL Lac's and other objects correlate strongly with the cosmic rays and will give quantitative 
predictions for the future experiments, in particular for the Pierre Auger Observatory. 
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